Quantum Stochastic Synchronization 
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We study within the spin-boson dynamics the synchronization of quantum tunneling with an 
external periodic driving signal. As a main result we find that at a sufficiently large system-bath 
coupling strength (Kondo parameter a > 1) the thermal noise plays a constructive role in yielding 
both a frequency and a phase synchronization in a symmetric two-level system. Such riveting syn- 
chronization occurs when the driving frequency supersedes the zero temperature tunneling rate. As 
an application evidencing the effect, we consider a charge transfer dynamics in molecular complexes. 
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The phenomenon of synchronization of nonlinear clas- 
sical systems with external driving signals, possibly even 
in the presence of randomness, has increasingly been 
gaining importance and growing interest over the last 
decade 0, 0, 0, 0- A particular intriguing example is 
crowd synchrony when pedestrians fall into steps with 
the intrinsic vibrations of a footbridge 0. Physically 
related to the synchronization phenomenon is the phe- 
nomenon of Stochastic Resonance |6| where the pres- 
ence of noise can manifestly boost the transduction of 
an information-carrying, weak signal. This effect is en- 
during ongoing vitality in diverse fields that span from 
physics to the life sciences. This fascinating phenomenon 
has notably been generalized into the quantum regime 
only recently 0,0111; where its experimental realization 
on the level of a nanomechanical quantum memory ele- 
ment is very feasible [ljj ■ The extension of noise-induced 
synchronization into the realm of quantum physics, how- 
ever, has not been considered thus far. This latter task 
presents a challenge which, apart from basic academic 
interest, also comprises great potential for nanoscience 
with promising applications ranging from quantum con- 
trol to quantum information processing. 

With this work we undertake a first step in this direc- 
tion. We consider the paradigm of a driven, symmetric 
quantum two-level system that is coupled to an Ohmic 
thermal bath with an exponential cutoff 0, [nj . When 
operating within the overdamped regime, we successfully 
reveal the existence of a noise-assisted quantum synchro- 
nization. There are two relevant energetic parameters 
characterizing the system-bath coupling: the reorganiza- 
tion energy A and the cutoff energy huj c which is associ- 
ated to the spectral width of the medium modes which 
couple to the tunneling particle. The two parameters 
can be combined into the dimcnsionlcss coupling strength 
a = X/(2Hw c ), i.e., the Kondo parameter. A third rel- 
evant energy scale constitutes the tunneling energy HA. 
The most relevant regime for our purposes corresponds 
to HA <C huj c , A and a > 1. 

Dissipative quantum tunneling changes radically the 
physics of synchronization. At zero temperature, the sys- 
tem can only tunnel towards its lowest energy state when 



a biasing dc-signal is applied. As the bias periodically 
changes its sign due to the action of a driving field, tun- 
neling makes the particle move periodically towards its 
corresponding lowest energy state, as long as the driving 
period is much longer than the typical tunneling time. 
Consequently, one expects that the system synchronizes 
with a driving periodic rectangular signal. By contrast, 
in the absence of classical thermal noise, synchronization 
in a overdamped bistable system driven by rectangular 
subthreshold signals does not occur, as the overbarrier 
transitions do not exist. 

Two interesting questions now emerge: What is the 
effect of quantum noise at finite temperature for syn- 
chronization? How does the time scale of the the exter- 
nal driving period impact this synchronization behavior? 
Quantum bath fluctuations at finite temperature surely 
will disturb the above mentioned perfect synchronization 
of the tunneling events as imposed by the rocking driving 
field. At the same time, when the driving period becomes 
shorter than the tunneling time, the zero temperature 
synchronization will also be weakened. Nonetheless, it 
might very well be possible that finite temperature quan- 
tum noise promotes and assists synchronization when the 
driving period decreases. We demonstrate below that 
this indeed is the case. 

To start, we consider the following archetype model of 
a dissipative, driven two-level system (TLS) 0, 0> 

m 

H{t) = ^e(t)a z + ^hAa x + ^a z Y,^(b]+b 1 ) 

j 
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Herein, the operators a z and a x denote the standard 
Pauli operators, e(t) is a time-dependent energy bias, 
and HA in Q is the tunneling matrix element. The 
boson operators b~^ and bj correspond to normal mode 
oscillators of the thermal bath with frequencies luj . The 
stochastic influence of the quantum thermal bath is cap- 
tured by an operator random force £(t) = ■ Kj(6^e l " j t + 
bje~ luJjt ). It can be characterized by the spectral den- 
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sity J(u>) — ('^/fi)J2j K 'j^( UJ ~ Uj)- We assume that 
J(u>) acquires the Ohmic form with an exponential cut- 
off, J(uj) = 2irhaujer ul / u]c . Due to the Gaussian statis- 
tics of a harmonic bath, the statistical properties of 
the quantum noise become completely determined by 
its equilibrium autocorrelation function (£(t)£(0))y = 

oc 

(h/ir) J J(oj)[coth(huj/2kBT) cos(wi) — i sin ut]duj. The 
o 

driven spin-boson Hamiltonian Q describes many sit- 
uations, such as, e.g., the electron transfer (ET) in a 
molecular dimer in azurin crystals |14| . Then the low- 
frequency molecular vibrations provide the bath and the 
time-dependent energy bias is given by e(i) = re£(t), 
where r is the tunneling distance, e denotes the charge 
transferred, and £(t) is the time-dependent, applied elec- 
tric field. 

In the incoherent tunneling regime, the populations 
of the localized states obey a nonstationary, Markovian 
dynamics. In the presence of a time-dependent bias, this 
description holds true for Ohmic friction at an arbitrary 
temperature if the tunneling coupling remains small, i.e. 
A <C u) c , and the coupling to the heat bath is sufficiently 
strong, a > 1/2 pjj. The populations p±(t) = (1 ± 
( ct «(*))t)/2 then obey the balance equations |3. UlL Il5| 



p p (t) = W-p(t)p-p{t) ~ W {t)pp{t) 



(2) 



with = ±, where the time-dependent relaxation rates 
Wp{t) are given by 

W±(t) = -A 2 J drexp[-Q'(r)] 

Q"{r) T \f t e{t')dt' 



X cos 
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within the Golden Rule approximation. The functions 
Q'(t) and Q"(t) in © denote the real and imaginary 
parts of 

Q(t) = — + ¥ j d tl j dt 2 (£(t 2 )£(0)) r , (4) 

wherein A = J °° cZwJ(aj)/(7ru;) is the bath reorganization 
energy 0. Here, A = 2ahu; c and the function Q(t) can 
be evaluated in closed analytical form 0, 0] , reading 



Q\t) = 2aln\ + 



T 2 (l + K ) 



-} (5) 



|r(i + K + ^ T t)| 2 . 

Q"{t) = 2a arctan(w c i) . (6) 

In Eq. (J5J, r(z) denotes the Gamma function, lot — 
ksT/H, and k = lot/lo c . Note that in the limit of an 
adiabatic driving varying on a time-scale T such that 
LD c T,aujTT 2> 1, the time-dependent transition rates 
W±(t) follow the instantaneous value of the bias e(t). 
In this limit, which is assumed throughout in the follow- 
ing, the relaxation rates W± (t) obey the Boltzmann rela- 
tion W+(t) = e e ^ /kBT W-(t). Furthermore, in the high- 
temperature limit ksT ^> fko c , Eq. ® reduces to a gen- 
eralized Marcus-Levich-Dogonadze form, i.e., W±(t) — 



(7r/2)SA 2 /V47rAfc B Texp[-(±e(t) - X) 2 /4Xk B T]. For 
ksT < fko c , explicit analytical expressions for the rates 
are generally not available, except for T = [13] ; typi- 
cally, those must be determined numerically from Eq. 

Using the Marcus-Levich-Dogonadze formula for an 
undriven molecular system, one can estimate the rel- 
evant parameter values. In particular, for ET occur- 
ring in azurin dimer the values are A = 0.25 eV and 
h A = 5 • 10 ~ 6 eV 0| . In molecular systems, the cutoff 
frequency of low-frequency molecular vibrations ranges 
between Hlo c = 5 — 20 meV. Choosing ho c — 12.5 meV 
yields a = 10 corresponding to a strongly incoherent, 
overdamped tunneling dynamics. 

It is worth noting that the present incoherent limit 
for the tunneling dynamics of a driven, dissipative TLS 
allows for an effective quasiclassical interpretation in 
terms of a classical random telegraph process which is 
inhomogeneous in time. Its transition rates, however, 
are governed by the manifestly quantum expressions de- 
tailed in Eq. ©. As such, our setup mimics the quan- 
tum analogue of a noisy, classical synchronization be- 
havior elaborated in Refs. [17], [l8|, [l9(. A thought- 
experimental setup is that of a particle tunneling be- 
tween two localized states at random times subject to 
an external periodic rectangular field with amplitude 
So and period T (frequency £1 — 2it/T). One counts 
the number of jumps N(t , to) within a time window 
[t ,t]. Following Refs. [H Hi 13, we introduce the 
random phase <j>(t,to) — wN(t,to), which increases by 
7r at each switching event (two subsequent switches cor- 
respond to a 27r-cycle with random duration) . We define 
the mean frequency and diffusion coefficients associated 
to the phase process as f2 p h := limt_ (00 (</>(i, to))/(t — to) 

and 2D ph := lim^ [(0 2 (Mo)> - (<£(Mo)) 2 ] /(* - *b), 
respectively. To evaluate these quantities, one can con- 
sider the joint probability Pp jTl (t) to be in the state (3 
at time t and to have n jumps within the time interval 
[to, t\. These probabilities can be obtained by integrating 
the multi-time probability densities of the corresponding 
stochastic trajectories [2(j. They satisfy the normaliza- 
tion condition Y^=o J2p=± Pp,n(t) — 1 an d obey a mas- 
ter equation resembling Eq. @: 

Pf3,n(t) = W-p{t)P-p, n -x{t) - W (t)P , n (t) (7) 

for n > 1 and P ,o(t) = -W [t)P 0jO {t) for n = 0. The 
populations of the states can be obtained as pp(t) = 
S^Lo Pfi,n(t) and the probability P n (t) to have n-jumps 
is then P n (t) = P+,„(i) + P-,n(t). Given P n (t), any 
moment (n k (t)) := J2^ =0 n k P n (t) (k = 1,2,...) can be 
obtained. Deriving explicit analytical expressions, how- 
ever, presents a nontrivial task. Fortunately, the first 
two moments are an exception to this rule. In particular, 
from Eqs. J7J and J3J, it follows that the phase frequency 
^ph(i) : = ir(d/dt){n(t)) can be expressed as 



fiph(t) = ir{W+(t) P+ (t) + W-(t)p-(t)] 



(8) 



This result coincides with one obtained earlier in Refs. 
[iH Il9| . For the averaged phase we have to)) = 
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fl p h(t')dt'. To obtain f2 p h, one must take the limit 
£l p h = lim^oo (4>(t, t Q ))/(t — t Q ). For a periodic driving 
the asymptotic behaviors for both pp(t) and S7 p i 1 (t) are 
periodic functions of time with period T. Thus, f2 p h = 
(l/T) £ dtSl§\t), where Sl$\t) is given by Eq. © 
with pp(t) replaced by the asymptotic solution Pp°^ (t) of 
the master equation (J2J), which formally is obtained by 
letting to — ► — oo. 

The calculation of the phase diffusion coefficient 
£> ph (t) := (n 2 /2)(d/dt) [(n 2 (t)) - (n(t)) 2 ] in closed form 
turns out to be intricate and lengthy, yielding 



2D ph (t) = ntt ph (t) 

xp 2 (t')exp 



2ir 2 5W(t) y p 



W{r)dT 



dt'W0) 



to 



(9) 



where 5W(t) = W+{t) - W-(t) and W(t) = W+(t) + 
W-{t). This remarkable exact result in Eq. © was 
originally derived in a different context in Ref. 0] 
within a slightly different approach. Note that D p h(t) 
still depends, like f2 p h(i), on the initial time tg. This de- 
pendence is asymptotically lost in the limit to — > — oo, 
where pp (t) is replaced by the asymptotic periodic solu- 
tion (t) . The evaluation of the mean diffusion coeffi- 
cient D p h proceeds similarly to f2 p h- For the considered 
case of rectangular driving, both, f2 p h and D p h are eval- 
uated to read: 



P h 



nW 



and 



2D 



1 - 5pL 



2vr 



4tanh(VFT/4) 



WT 



(10) 



P h 



7T O ph - _<5p 4 q tanh d (WT/4) 



-%f S Pl q (! - S P 2 J { 12 tanh ( WT / 4 ) 
-WT [l + 2sech 2 (WT/4)] }, 



(11) 



where 14 7 denotes the sum of the forward and backward 
rates in Eq. 10 for a fixed value of the field £q, i.e., 
for e(t) = eo = er£$, and (5p cq = tanh(eo/(2fcsT)) is 
the absolute value of the difference of the equilibrium 
populations. The inverse Fano factor of the counting 
process R :— 7rf7 p i 1 /(2Z3 p h) provides a reliable quality 
measure of synchronization 0, 0] ■ 

The quantum features of the so derived synchroniza- 
tion quantities are rooted in the quantum rate expres- 
sions entering Eqs. (|10f> and Hll|) . By analogy with the 
findings for quantum stochastic resonance (QSR) in sym- 
metric quantum systems one could expect that there 
is no thermal noise assisted synchronization for a < 1. 
Indeed, we could not identify noise-assisted synchroniza- 
tion in this parameter regime. In contrast, employing the 
reasoning put forward for QSR within a > 1 in Ref. [8|, 
one then supposes quantum synchronization to emerge in 
this latter regime. 



Indeed, for a > 1 the synchronization scenario 
depends sensitively on the value of the driving 
frequency. As discussed above, for driving time 
scales with a frequency much smaller than the 
tunneling rate at T = 0, i.e., ft <C Wt=o — 
nA 2 [2uj c T(2a)}- 1 [eo/(hij c )} 2a - 1 exp[-eo/(hLj c )}, O ph 
matches the external frequency for sufficiently low 
temperatures. This is depicted in the upper left panel 
of Fig. ^ where in addition the behavior of -D p h is 
presented. Notice that quantum synchronization at 
sufficiently low temperature indeed is impressively well 
accomplished, as it is reflected by the large values of 
the R- factor (see the lower right panel of Fig. As 
the temperature increases above a certain critical-like 
value, however, the R- factor diminishes and the quality 
of synchronization deteriorates. 

With increasing driving frequency such that the con- 
dition Q -C Wt=o is no longer being met, we enter the 
regime of thermally induced phase synchronization in the 
quantum regime, as depicted in the upper right and lower 
left panels in Fig. ^ We observe a range of temperature 
where the external frequency and f2 p h coincide. More- 
over, the diffusion coefficient exhibits a pronounced cusp- 
like minimum. The results plotted in those two panels 
convincingly illustrate the constructive role of quantum 
thermal noise for quantum synchronization. Notice that 
with increasing driving frequency the quality factor as- 
sociated to the synchronization effect starts to diminish, 
see the lower right panel in Fig.^ Moreover, the temper- 
ature range where frequency locking is observed shrinks. 
The quality factor R displays a cusp-like feature similar 
to the one discussed recently by Park and Lai [2l| within 
the context of noisy, classical synchronization. Upon fur- 
ther increasing the driving frequency, the tunneling rate 
is too slow compared to the driving frequency for the tun- 
neling dynamics to follow the external oscillations. Thus, 
quantum synchronization is lost (not depicted). 

In conclusion, we have discovered the existence of 
a quantum stochastic synchronization in an externally 
driven spin-boson system which undergoes tunneling 
transitions between two states. We exemplified the phe- 
nomenon for non-adiabatic charge transfer in molecular 
complexes. The estimates of parameter values for an ex- 
perimental test of our theoretical predictions are: driving 
frequencies O ~ 10 _1 — 10 3 s _1 , electrical field strength 
£ - 10 4 - 10 5 V/cm, and temperatures T ~ 20 - 100 
K, which are readily achieved in the laboratory. The 
main quantum features of the discussed synchronization 
phenomenon are robust and they are not critically depen- 
dent on the details of the underlying dissipation mech- 
anism. This is so because the phenomenon of quantum 
synchronization is primarily based on the existence of a 
low-temperature limit of the tunneling rates. We, con- 
sequently, expect our very general results to be useful in 
other contexts such as for optimizing the noisy dynamics 
of nano- mechanical systems in the quantum regime |10| . 
The authors are confident that these riveting findings for 
quantum synchronization will spur experimental interests 
for diverse systems that involve a controllable quantum 
tunneling between distinct quantum states. 




FIG. 1: Mean jump frequency fi p h, phase diffusion coefficient D p h (upper right, upper left and bottom left panels) and the 
synchronization quality factor R (bottom right panel) versus the scaled temperature n — kBT/hui c for three values of the 
driving frequency fl. The parameter values used are: a = f 0, A = 4 ■ 10 -4 , eo = 5. Frequencies are scaled with lo c and energies 
with huj c - For huj c = 12.5 meV, ui c ~ 1.9 • 10 13 1/sec, and r = 14.9 A (taken to match the charge transfer in molecules from 
Ref. |3) yields £o ~ 4.19 ■ 10 4 V/cm. The value « = 1 corresponds approximately to 145 K. 
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